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PHOTOELASTIC ANALOGY FOR NON-HOMOGENEOUS FOUNDATIONS 


Allen Curtis* 
and 
Richart, Jr.**, Assoc. ASCE 


Synopsis 


This paper deals with the stress distributions semi-infinite body 
which ratio constant, but the modulus elasticity varies 
with depth. The general equations elasticity for such point-isotropic 
material are derived and the validity and limitations published solutions 
are discussed. There shown analogy between the stresses 
constrained slice, which has modulus elasticity variable with depth, 
and the stresses plate whose thickness varies the same manner, 
providing the state stress the plate may considered two- 
dimensional. Thus the stress distributions beneath infinite line strip 
loads such foundation may determined experimentally the use 
photoelasticity. The photoelastic analogy confirmed experiments 
and its use discussed detail. 


Introduction 


For certain problems the field soil mechanics, the stress-strain 
relations the mass soil agree reasonably well with Hocke’s law. For 
such conditions, the theory elasticity will give satisfactory results for 
stresses and displacements the properties the theoretical material 
are similar the physical properties the soil mass.! the soil mass 
homogeneous and isotropic well elastic, solution involving only 
two elastic will applicable. Such solutions that 
for concentrated normal load, and that Cerruti’, for 
concentrated tangential load, each acting the free surface semi- 
infinite body, have often been used satisfactorily for the study stresses 
and deformations beneath foundations. suitable integration the 
equations for stresses, given these fundamental problems, has been 
possible evaluate the stress distributions beneath line, strip, and areas 
the surface which are loaded. 

However, the materials encountered soil mechanics problems are 
seldom both homogeneous and Consequently, actual results are 
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often variance with the predictions from theory. result this 
discrepancy, numerous studies have been made, still based the theory 

elasticity, but for which the idealized elastic material assumed 
have properties more nearly like those soils. For various analyses, 

the elastic material was assumed be: (a)homogeneous but aelotropic, 

that is, having different elastic properties the horizontal and vertical 
directions, (b) made layered systems, each layer being homogeneous 
and isotropic, (c) point isotropic, which the material assumed 

isotropic each point, but non-homogeneous since the elastic con- 
stants vary with their location. 

This investigation concerned with the analyses based assumption 
point-isotropic materials for which the Poisson’s ratio constant 
throughout, but for which the modulus elasticity varies with depth. First 
the general equations elasticity are derived for this particular type 
point-isotropic and these are used evaluate the existing 
theoretical studies the problem. Then comparing the general equa- 
tions elasticity with those equations which represent the stress dis- 
tribution plate variable thickness, analogy established which 
permits the analysis stress distributions the point-isotropic material 
the use photoelastic model which has variable thickness. 

the purpose this paper establish the validity the photo- 
elastic analogy, and present study its usefulness and limitations 
aid for the study stress distributions point-isotropic materials. 


Notation 

The letter symbols used this paper generally conform those pre- 
scribed the Manual, ‘‘Soil Mechanics Nomenclature’’, published the 
American Society Civil Engineers, 1941. 

Also, the interest brevity, the term will used 
the following pages indicate solution for the stress distribution 
semi-infinite body, due surface loading. 

the following discussions, compressive stresses are considered 
positive customary the study soil mechanics. Figure shows 
the coordinate systems used and Figures 2(a) 2(e) show the directions 
positive stresses the faces infinitesimal element material, 
defined these coordinate systems. 


Equations Elasticity for Point-Isotropic Materials 


The theory elasticity Hooke’s law, which may expressed 
general form the ‘‘Each the six components stress 
any point body linear function the six components strain 
the point.’’ This law then gives possible total independent 
elastic constants each point the material. However, the material 
assumed isotropic, that is, the elastic properties are independent 
rotation coordinate axes, only two independent elastic constants 
are obtained any point. If, addition, the assumption made that the 
material homogeneous, only two elastic constants are obtained through- 
out the body. Since the mathematical procedures are greatly simplified 
using only two elastic constants, the greater part the work the 


University Press, London, 1934, 97. 
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theory elasticity has been done the assumptions that the material 
isotropic and homogeneous. 

For the point-isotropic material, only two independent elastic constants 
are obtained each point the body. However, these two quantities may 
each have different values different locations within the elastic body. 
The values Poisson’s ratio, and the modulus elasticity, will 
used this study the two independent elastic constants. will 
assumed that Poisson’s ratio has the same value throughout the body, but 
that the modulus elasticity function depth only. Since the modulus 
term may used equally well the subsequent equations. 

These foregoing assumptions regarding the elastic constants way 
affect the usual stress-strain relations, since these are concerned with 
the conditions point the body. Also the strain-displacement equa- 
tions are unaffected the assumed variation the elastic constants. 

The equilibrium equations terms stress express the equations 
statics which must satisfied throughout body. Thus again, variations 
elastic constants cannot affect the equations this form. terms 
stresses, the equations equilibrium may written, 


for Cartesian coordinates. The (x, the right the equation indi- 
cate that similar equations, for the other coordinate directions, may 
written cyclic interchange subscripts. 

the type problem under consideration, the only body forces en- 
countered will the gravity loading the soil. However, will shown 
that the differential equations satisfied this point-isotropic medium 
are linear, that the law superposition valid. Thus, the effects due 
gravity loading may added separately any solution obtained; con- 
sequently the body forces will neglected the following pages. 

addition the equilibrium equations, boundary conditions, and stress- 
strain relations complete solution problem use the theory 
elasticity must also satisfy the equations compatibility. These equations, 
when satisfied, guarantee that the displacements the body are single- 
valued every point. terms strains, the compatibility equations for 
three-dimensional system are, 


im 
a 


These equations are independent the elastic constants, and 
this form. However, one convenient method attack for solving prob- 
lems elasticity transform Eqs. (2) and (3) into terms stresses, 
making use the stress-strain relations and equilibrium 
When the boundary conditions are known terms stresses, the 
case for the type problem under consideration, this seems pre- 
ferable. 

accordance with the usual practice, the z-direction has been chosen 
that is, 


(4) 


Substituting the stress-strain relations into Eqs. (2) and (3), and using 
the equilibrium equations eliminate shearing stresses, the three- 
dimensional equations compatibility terms stresses become: 


(5) 


Hill 1951, 230. 
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operator, denotes the sum the principal stresses, and designates 
the natural logarithm the shear modulus. 
Adding Eqs. (5), (6), and (7), and again using the equilibrium equations 
eliminate shear stresses, the following equation obtained: 


(8) 
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evident that constant, the above compatibility equations 
reduce those given Timoshenko and Moreover, all these 
equations are linear that the stresses and displacements due different 
boundary conditions may superposed. 

the system which corresponds plane strain the plane, the 
compatibility relations are reduced from six one, since 


substituting these relations into Eqs. (2) and (3), and using the stress- 
Strain and equilibrium equations, the compatibility equation terms 
stress determined follows for plane strain: 


where 


For the case plane stress, slice considered state 
plane stress the faces the slice are free surface tractions, and 
thin enough comparison with its other dimensions that vari- 
ations stress through the thickness may considered negligible. 

This condition requires that for plane stress the plane, 


Under these assumptions, the compatibility relations reduce to, 


DEx = ote 


See footnote 231. 
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Eqs. (16) are satisfied 


customary solve plane stress problems neglecting Eqs. (16) 
which also involves disregarding two shear stress-strain relations which re- 
late Eqs. (14) and (16). has shown that the error due this 
simplification small the change the normal stresses moderate 
distance equal the thickness the slice. addition, integrating 
the correction the stresses obtained this manner, across the thickness 
the plate, the result zero. Consequently photoelastic work the 
fringe patterns will agree with this average stress distribution. 

The compatibility equation for the case plane stress the 
plane may transformed into terms stress give: 


The differential equations which have been derived above for point- 
isotropic material could put many forms, for example terms 
Airy’s stress function for the two-dimensional cases. evident that 
the variable modulus elasticity greatly complicates the equations 
satisfied, and there seems little hope for mathematical solutions 
practical value except for special forms some which will con- 
sidered the following section. 

Inspection these equations shows the following points. First, all the 
equations are linear, that the solutions for different surface loadings 
displacements may superposed. Thus the solutions for concen- 
trated force line load infinite length are again fundamental solutions 
for the three- and two-dimensional problems respectively. Secondly, 
the rate change (or which influences the stress distri- 
bution, rather than the actual value rate change (or E). Finally, 
the solutions two-dimensional problems depend the value Poisson’s 
ratio, contrast the solutions ordinarily obtained under the assumption 
constant modulus. 


“On the Calculation Three-Dimensional Corrections for the Two- 
Dimensional Theory Plane Stress,’’ Reissner, Proc. 15th 
Semi-Annual Eastern Photoelasticity Conference, June 1942, 23. 


z 
(18) 
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Previous Solutions for Point-Isotropic Materials 


Since the solutions for homogeneous isotropic materials are well known, 
they will serve convenient standards comparison for solutions based 
different types materials. For normal point load acting the 
surface the semi-infinite homogeneous isotropic body, 
equations for the stresses are available various books dealing with the 
theory 

Boussinesq’s equations for stresses, for the condition that 1/2, 
correspond what has been called ‘‘rectilinear’’ transmission stress. 
this meant that the direction any ray emanating from the point 
load application the direction the only non-vanishing principal 
stress. However, the term ‘‘radial transmission stress’’ seems 
more descriptive and will used the following discussion designate 
this particular stress distribution. For this stress distribution, the only 
principal stress not equal zero, is: 


interesting note that for any value Poisson’s ratio the result- 
ant the normal and shear stresses element horizontal plane, 
passes through the point application the load. This evident from 
the relation, 


(20) 


has shown that for any value Poisson’s ratio, the direction 
the maximum principal stress does not vary appreciably from the radial 
(R) direction until the angle becomes large. this region however, the 
stresses are now small. Accordingly assumed radial transmission 
stress approximation for any value Poisson’s ratio, and for any 

The equation for the equilibrium the region beneath the load estab- 
lishes stress distribution, 


1934. 
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where arbitrary parameter. Since this the only principal stress 
not equal zero, the stresses expressed terms cylindrical coordin- 
ates are obtained directly as, 


Oz = 


For values n>2, the above stresses satisfy the differential equations 
equilibrium and the boundary conditions. 

obtain the appropriate value Frohlich determined the value 
which made the strain energy minimum throughout hemispherical shell, 
radii and (R, R,), beneath the load. introducing the expres- 
sion for the variation the modulus, 


(23) 


the strain energy the region was found minimum for 
For constant modulus, and substitution this into Eq. (21) gives 
the stress distribution for 1/2, (Eq. (19)). For modulus 
increasing linearly with depth, 

The value obtained above independent the inner and outer radii 
the shell; consequently the result valid for semi-infinite body, pro- 
viding small region about the singularity under the load excluded 
the usual way. 

Frohlich has thus obtained stress distribution which satisfies the 
boundary conditions, the differential equations equilibrium, and makes 
the strain energy minimum for the assumed form the stress distrib- 
ution. However, this does not mean that the solution correct, since 
the compatibility relations have not been fulfilled. The minimizing the 
strain energy does not replace these relations; effect, makes the 
deviation from the correct solution small possible for the assumed 
stress distribution. 

The form Eq. (21) has suggested the interpretation the quantity, 

‘‘stress concentration factor.’’ This concept permits evaluation 
the influence non-homogeneity the material upon the stress distri- 
bution beneath the load. has integrated Eq. (21), obtain the 
values the axis symmetry due uniform and parabolic 
intensity pressure distribute? over circular area the surface. 
used values equal and and showed that the stresses computed 
using and the parabolic surface load distribution agreed best with 
experimental data, depths greater than the radius the loading area. 
However, did not mention the elastic constants the material used, and 
concluded that the equations elasticity must modified before they can 


Transactions, ASCE, 101, 1936, 1072. 
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applied soils. will shown later, the elastic constants the 
material determine the value which must used Frohlich’s equa- 
tions give correct solution. 

the discussions Cummings’ paper, various investigations presented 
correlations between this quantity, and properties soils. Krynine 
derived expression, 


(24) 


2 


which the coefficient earth pressure rest, defined the ratio 
horizontal vertical pressure point. This ratio may also deter- 
mined as, 


ta se 25 
(25) 


any point semi-infinite elastic, isotropic mass material which 
loaded its own weight. Thus, Eq. (24) can rewritten as, 


Krynine, Newmark, and Casagrande each pointed out the possibil- 
ity that might vary with depth actual soils. Also was noted that layers 
different stiffness affected this stress concentration factor. Thus the 
modulus elasticity the material, and more specifically its mode 
variation with depth, well Poisson’s ratio, were believed have 
significant correlation with the factor, 

Further studies stresses non-homogeneous foundations were made 
Ohde. examined Frohlich’s solution with respect the question 
the compatibility relations, and employed equation compatibility cor- 
responding Eq. (11). considered moduli elasticity given by: 


(27) 


For these moduli, Eq. (11) becomes, respectively, 


(29) 


Theorie der Druckverteilung Baugrund,’’ Ohde, Der 
Bauingenieur, Vol. 20, 1939, 451. 
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which Ohde stated that Eq. (30) satisfied Frohlich’s solu- 


tion (Eq. (21)) 


For correct solution, all six compatibility relations must satisfied. 
Substitution the stresses, obtained from Frohlich’s solution, into Eqs. 
(5) through (10) shows that they are satisfied for modulus given Eq. 
(28), Eq. (31) satisfied. 

The significant point concerning Eq. (31) that defines definite 
relationship between ‘‘stress concentration factor,’’ (n), 
Poisson’s ratio, and the exponent which describes the manner vari- 
ation the elastic modulus with depth. Thus for constant modulus, 

solution identical with that Boussinesq. 
modulus increasing linearly with depth defined, and correct results are 
obtained from solution only Poisson’s ratio 1/3. 

Frohlich has also obtained solutions (See Table similar manner 
for the two-dimensional problem. Calculations show that this solution 
again satisfies the compatibility relation for plane strain (Eq. (13)) under 
the same conditions for the case. This ex- 
pected since the two-dimensional solution merely integration the 
three-dimensional solution. further integrations, solutions correspond- 
ing any given normal surface loading can obtained. 

Ohde also derived the correct solution for concentrated force, 
applied the direction the x-axis point the surface the semi- 
infinite body. The stress distribution defined by: 


n-3 


For correct solution, the value used Eq. (32) must determined 
from Eq. (31) the modulus elasticity defined Eq. (28). 
the above equation becomes identical the one given Cerruti. 

The solution derived Frohlich for the condition involving modulus 
elasticity increasing linearly with depth kz), has been confirmed 
Borowicka.!2 assuming the stresses and displacements terms 
infinite power series with unknown constant coefficients, obtained solu- 
tion for the case the point load which may evaluated for any value 


12. “Die Druckausbreitung Halbraum bei linear zunehmenden Elastizi- 
tatsmodul” Borowicka, Ing. Archiv., Vol. 14, 1943, p.75. 
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Poisson’s ratio. For two values the series degenerate single 
terms. For Frohlich’s solution obtained, and for 1/2, the 
result identical with the Boussinesq solution. Since Borowicka worked 
directly with displacements, that is, the basic equation elasticity, 
was unnecessary use the compatibility equations. Substitution the 
stresses obtained using 1/3 1/2 into the compatibility equa- 
tions shows that these are satisfied his solution. Borowicka’s results, 
shown Table indicate that the values the vertical stress the axis 
loading, depend very much upon the value Poisson’s ratio. 
change Poisson’s ratio from 1/2 zero, corresponding increase 
stress about 76% occurs, any given depth. Consequently, the use 
solution may give considerable error the conditions Eq. 
(31) are not satisfied. 

further study, Borowicka presented solution for the two- 
dimensional problem for which the modulus elasticity the material 


was assumed be, 


(33) 


which corresponds infinitely rigid top surface. For this case, the 
last term the compatibility relation for plane strain (Eq. 13) equal 
zero, that the stress distribution independent Poisson’s ratio. 
The stress distribution radial, and the only principal stress not equal 


zero given by, 


(34) 


Along the axis loading, this gives the same value the 
solution for constant 

The condition that Eq. (31) must satisfied the value used 
Frohlich’s solution, does not imply that must integer. However, 
must equal to, greater than since Poisson’s ratio cannot 
greater than 1/2. The expression for the modulus elasticity, 


(35) 


requires the use 3.5, and 0.4 the Frohlich equation order 
arrive correct solution. The modulus described Eq. (35) has 

been used approximation for the study stress distributions beneath 
footing granular soils. solution for stresses and displacements, 
based this variation the elastic modulus, was given Hruban. 


13. Druckausbreitung einer Halbscheibe bei mit der Tiefe abnehmen- 
den Elastizstatsmodul,’’ Borowicka, Oest. Ing. Archiv., Vol. 
1948, 360. 

14. See footnote for discussion this solution. 

15. Influence the Heterogeneity the Soil Mass its Deformation’’, 
Hruban, Proceedings the Second International Conference 

Soil Mechanics and Foundation Engineering, Rotterdam, 1948, Vol. 


123. 
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considered both normal and tangential loadings, acting point, strip, 
Circular area the surface the body. 

The preceding section has established the validity and/or limitations 
the known solutions for point isotropic materials. addition, has 
provided some check the derivations the general equations elas- 
ticity, inasmuch Ohde and Borowicka arrived the same equations 
for particular types variable moduli. 

The available results for the vertical stress, are summarized 
Tables and 


The Photoelastic Analogy 


the previous sections, the effects the stress distributions due 
modulus elasticity varying with depth have been examined from 
mathematical point view, and the differential equations which must 
have been derived. 

instructive consider the problem from physical viewpoint 
also. Data from soil tests indicate that elastic soils often exhibit modu- 
lus which increases with depth. Such soils may perhaps considered 
made infinite number thin horizontal layers homogeneous 
isotropic materials. With increasing depth, each layer somewhat more 
rigid that given uniform pressure two adjacent layers will pro- 
duce smaller contraction the lower layer. Such physical concept 


methods estimating the stress distribution non-homogeneous founda- 
tions. these methods the vertical thicknesses layers were adjusted 
proportion their relative stiffnesses, and then the Boussinesq solu- 


tion was applied this transformed material. 

similar physical interpretation the problem may gained con- 
sidering thin slice from semi-infinite body homogeneous, isotropic 
material which loaded along region the free surface. the slice 
has constant thickness, then the stress distribution will correspond 
that produced infinite line strip loading the semi-infinite body. 
Now let the thickness the slice vary function depth only, and 
assume increases with depth. the thickness changes smoothly, and 
not too rapidly, then reasonable expect that the stress distribution 
may considered two-dimensional, and that the variations stress 
through the thickness are negligible. The slice may thus considered 
infinite number thin strips, each strip being somewhat wider with 
increasing depth. uniform pressure applied over the top surface 
such the value stress and accompanying vertical contractions 
will decrease with depth. Thus there appears physical analogy 
between the soil with varying modulus, and the slice with varying thick- 
ness. will now Shown that this relationship holds 
thus establishing method for determining the stress distribution 
point-isotropic material the use photoelasticity. 

The stress distribution semi-infinite body produced line 
strip load the surface analyzed problem plane strain. How- 
ever, photoelastic model represents problem plane stress. order 
use the results from photoelastic tests describe stress distributions 


und Bettungsziffer’’, Terzaghi, Oesterreichische 
Bauzeitung, No. 25, 1932. 
17. Footnote, 10, see Discussion Newmark, 1115. 
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for case plane strain, necessary use the plane stress plane 
Strain, free slice constrained slice, analogy which has been discussed 
Westergaard Mindlin and Salvadori. 

For homogeneous isotropic materials, this analogy states that any prob- 
lem plane stress may considered problem plane strain, 


vice versa, provided that suitable change the elastic constants made. 


the region under consideration, (a) not subjected body forces, 
subjected body forces which are harmonic, (b) simply connected, and 
(c) not subjected thermal stresses, then the stress distribution 
independent the elastic constants, and the same for free and con- 
strained slices. The problem under consideration using point-isotropic 
materials satisfies conditions (a), (b), and (c) above, but not independ- 
ent the elastic constants due the change the compatibility equation. 
(see Eqs. (13) and (18)). therefore necessary determine what 
changes elastic constants are necessary for the analogy hold for 
point-isotropic materials. 

Consider plane stress system with elastic constants E', 
plane strain system with constants and Then may 
the substitution 


and 
Shown that 


and 
(36) 


into the stress-strain relations for plane stress yields the equations for 
plane strain. Similarly, the substitution 


(37) 


into the stress-strain relations for plane strain yields the equations for 
plane stress. The equilibrium equations terms stress are identically 
equal and inspection Eqs. (13) and (18), seen that the compatibil- 
ity equations also become identical the above substitutions are made. 
Thus free slice and constrained slice, which have the same physical 
dimensions and the same boundary loading, will have the same stress dis- 
tribution the Poisson’s ratios the two materials are related shown 
above and the moduli elasticity are directly proportional. 
interest note that the free slice constrained slice relations between the 


18. ‘‘Graphostatics Stress Functions’’, Westergaard, Trans- 


actions ASME, Vol. 56, No. 1934, 141. 

19. Experimental Stress Analysis’’, edited Hetényi, 
Chapter VII ‘‘Analogies’’, Mindlin and Salvadori, 
Wiley and Sons, 1950. 
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elastic constants for this point-isotropic material are the same those 
for the case homogeneous isotropic materials which the stress dis- 
tribution dependent the elastic constants. (see footnote 17, 762). 
From Eq. (36) seen that practical use this free slice-constrained 
slice analogy limited certain range values for Poisson’s ratio. 
the value the plane strain system 1/3, then the corresponding value 
for plane stress 1/2. Since impossible for material have 
value Poisson’s ratio greater than 1/2, the limiting value for the 
analogy 1/3, (plane strain), corresponding 1/2, (plane stress). 
Now consider semi-infinite thin free slice, subjected given bound- 
ary loading, whose thickness (h) function the depth (z). The thick- 
ness assumed vary gradually, that the state stress may con- 
sidered two-dimensional. The equations equilibrium infinitesimal 
element this slice may written as, 


and 2 z 


The stress-strain relations may written, 


and the compatibility relation given Eq. (15). use Eqs. (25), 
(38), and (39), the compatibility equation for the variable thickness slice 
given terms stresses as: 


comparing Eq. (40) with Eqs. (13) and (18), and comparing the 
equilibrium equations, seen that the two systems are mathematically 
identical, except for the boundary conditions, one identifies: 


(40) 


age 
(38) 
(29) 
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Therefore, semi-infinite free slice with varying thickness loaded 
the top surface, the stress multiplied the thickness given point 
equal the stress corresponding point semi-infinite constrained 
slice varying modulus, loaded similarly, providing that the relations 
Eq. (41) are maintained. Thus the analogy proved mathematically. 
should noted that the compatibility relation, merely the deriva- 
tives the and terms which must correspond, that the thick- 
ness the free slice may scaled geometrically any desired value. 

Since obviously impossible test semi-infinite slice, either 
photoelastically otherwise, finite model must used. and 
Frocht2! have used relatively small sheets Bakelite photoelastic 
models check the Flamant solution for the semi-infinite plate. Away 
from the edges the model, the results checked very well with theory. 
Thus, comparing the stress distribution two plates the same size 
under the same load, one plate having constant thickness and the other 
having variable thickness, appears possible determine use 
this analogy the effect variable modulus has upon the stress distribution 
caused infinite line strip loads. 

The photoelastic method provides means determining each point 
the model, (a) the difference the principal stresses, from the fringe 
patterns, and (b) the direction the principal stresses, from the iso- 
clinics. For complete study the stress distribution, also necessary 
determine the sum the principal stresses each point. This may 
done any one the usual methods. 

For this type problem appears that numerical method will prove 
the most expedient since the models will have rectangular boundaries. 
The boundary conditions will consist load applied point, distri- 
buted over small area, the top edge the plate and known as- 
sumed distribution stress the bottom edge. The vertical edges 
the plate will free from external loads. comparison made be- 
tween the stress distribution the plate constant thickness and that 
the plate whose thickness varies with depth. 

For the case the constant thickness plate, the photoelastic fringe 
pattern will determine the values the difference the principal stresses 
each point the model the relation: 


20. ‘‘Distribution Stresses Partial Compression and...’’, 
Durelli, Proc. 13th Semi-annual Eastern Photoelasticity Conference, 
June 1941, 15. 

21. ‘‘Photoelasticity’’, Frocht, Vol. II, Wiley and Sons, Inc., 


Y., 1948. 
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the material fringe constant, 

the thickness the model. 
the free boundary, equals that the fringe pattern also 
determines the boundary stress values. The singularity due the concen- 
trated load the top surface may eliminated the usual way con- 
sidering small semi-circular groove cut out under the load. Then 
the surface the groove assumed subjected equipollent 
radial stress distribution given the Flamant solution. The stress 
distribution along the base the model will determined the method 
supporting the model. The boundary values established 
above, and the compatibility relation, 


(43) 


which results from Eq. (40) for constant thickness plate, are sufficient 
for the solution the stress distribution the iteration relaxation 
methods. The finite difference net for use with such numerical methods 
shown Fig. 

Now consider the case the photoelastic model having variable thick- 
ness. Eq. (42) may rewritten as, 


and the compatability equation satisfied given Eq. (40) 
substituted for The boundary values h(p+q) are established 

the same manner for the constant thickness model. order solve 
Eq. (40), necessary express interms h(p+q) and other 
known quantities. introducing the angle, between the z-axis and the 
direction the principal stress, Eq. (40) becomes, 


(45) 


The only unknown Eq. (45) h(p+q) interior points, since the values 
are known from the isoclinics, and h(p-q) known from the fringe 
pattern. Thus, Eq. (45) may solved numerically and the complete stress 
distribution the model may obtained. 
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Figure shows Eq. (45) the finite difference form for square net. 
the variation thickness not easily expressed analytically, the re- 
quired derivatives may also obtained graphically the use finite 
differences. 


Selection Photoelastic Model 


The analogy has been proved mathematically, but still necessary 
evaluate its suitability tool for experimental stress studies. 

The first consideration must the choice photoelastic model which 
corresponds analogy given variation elastic modulus the semi- 
infinite body. the known solutions discussed, the only one which valid 
for any given value Poisson’s ratio Borowicka’s solution using kz. 
The photoelastic model which represents this modulus has thickness which 
corresponds kz, wedge. The obvious difficulties encountered 
applying load the sharp edge wedge are partially offset the simple 
machining operations involved preparation the model. 

After consideration the photoelastic behavior, machinability, and 
availability the various materials ordinarily used for photoelastic models, 
Bakelite 61-893 was selected for the wedge. The material fringe value 
was found use standard tension specimens 44.3 psi. shear per 
fringe per inch thickness for the 1/2 inch thick plate, and 41.5 psi. shear 
per fringe per inch for the 1/4 inch thick plate. Poisson’s ratio this 
material was found 0.366 from lateral and longitudinal strain measure- 
ments strip. The longitudinal strains were measured with two 
inch gage length Huggenberger strain and the lateral strains were 
measured use Vose Interferomete 2strain gage. 

After selection suitable material for the photoelastic model, the 
physical dimensions the model were chosen. Since the test was simu- 
late the behavior semi-infinite plate, obvious that the height-width 
ratio finite model will have important effect the stress distribu- 
tion. Previous have shown that the theoretical fringe pattern pro- 
duced line load acting semi-infinite plate, consists circles tan- 
gent the top surface the point loading. Thus theoretical fringe 
would tangent the four sides square cut out the plate. However, 
photoelastic model usually supported along the base, which could con- 
sidered act extending the model. addition, the vertical sides 
loading can expected have greater effect the stresses the neigh- 
borhood the load than does proximity the base. This leads model 
wider than high. series tests was conducted determine the most 
desirable height-width ratio for plate loaded concentrated load and 
supported distributed reaction. From these tests height-width ratio 
0.60 was chosen. 

order derive full information from use the photoelastic anal- 
ogy, very important that the ‘‘stress effects due the 
variable thickness are separated from those effects introduced the bound- 
ary conditions geometry and stress distribution. analysis the 
effect these boundary conditions was made during this investigation. 

The next step determine suitable height and thickness for the 


22. application the Interferometer Strain Gage 
Vose, Journal Applied Mechanics, Trans. ASME, Vol. 
No. September 1935, A-99. 

23. See Footnote 85. 
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model. The overall size usually limited the size the polariscope 
field, unless movable loading frame devised that overlapping 
photograph technique could used. However, the primary consideration 
design model that the best possible fringe count obtained over 

the whole area. Increasing the thickness increases the fringe count, for 
given stress, but this simple answer counteracted the necessity 
maintaining thin enough model preserve two-dimensional stress sys- 
tem. wedge angle about 14° was found satisfactory for such 
models. The useful height the wedge determined the fringe count 
which may obtained. From Eqs. (41) and (44), seen that the maxi- 
mum shear for the analog given the value (nf) each point the 
model. Thus the fringes, lines constant are expected approxi- 
mately circular for the plate, and approximately the same size for 
the same total load. The total load the wedge analogous the line 
load the point-isotropic medium. However, since the edge the wedge 
can support only rather small load, the height the model limited 

value such that fringes are obtained over the whole model. Preliminary 
experimentation lead the choice height 0.90 inches. 

With regard the top surface the model, also obvious that 
thickness zero would develop infinite stress under any load. Consequent- 
ly, the edge would chip out, flow plastically, even under small loads. 
Also this chipping effect would likely occur during machining the Bake- 
lite doubtful such edge could machined, even were 
desired. Preliminary experiments showed that top thickness about 
0.020 inches was the best compromise between desire and possibility. 


Experimental Procedure 


The wedge was loaded along its plane symmetry concentrated 
force acting the top edge. The lower edge the model rested several 
thicknesses cardboard. The cardboard rested Plexiglas block which 
turn was supported steel block. Steel rollers were placed between 
the steel block and the steel base plate eliminate any horizontal 
force which might introduced the loading device. The method 
loading the models shown diagrammatically Fig. and the actual test 
set-up shown Fig. The load was applied through lever system and 
was indicated the deformation series iso-elastic springs. This 
loading frame was calibrated both tension and compression before start- 
ing the tests. 

Since the wedge was loaded along its axis symmetry, its faces were 
not perpendicular the optical axis the polariscope. Thus was neces- 
sary immerse the wedge liquid bath which had the same index 
refraction that the model material. satisfactory liquid was prepared 
mixing one part Bayoi White oil with two parts Halowax 1007 oil. 
The container for the oil was made from sheets 1/8 inch and 1/2 inch 
Plexiglas and shown Fig. 

order increase the fringe count obtainable for given load per 
wedge, two wedges were placed parallel and loaded single pin. This 
procedure involved some question the distribution load each 
wedge, although the results obtained were found agree well with later 
tests single wedges. However, the difficulties involved equalizing 
the loads, and also those alignment made this set-up impractical; con- 
sequently single wedge was used for all remaining tests. 
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The thin top edge the wedge caused large stresses appear this 
region, even when relatively light loads acted. order obtain fringe 
values close the load before the pin dug into the edge, photographs were 
taken the fringe pattern for series increasing loads. Thus the values 
obtained are most accurate the top portion the model for small values 
load and more accurate lower depths for large values load. 

Tests were also run plates constant thickness, evaluate the 
method support the base. plate Bakelite 1.80 inches high 1.50 
inches wide was loaded concentrated loads opposite edges, thus simu- 
lating the problem layer finite thickness slippery rigid base, 
(Fig. 6c). the next test, block was loaded shown Fig. 6d. This 
test analogous the problem two elastic layers rigid base, assum- 
ing (a) the interface between elastic layers perfectly adhesive, and (b) the 
surface the rigid base slippery, thus developing shear stress. 
the increase thickness, the stress distribution certainly not two-dimen- 
sional, but small region each side the jump thickness ex- 
cluded, the stress distribution analogous the layered system problem. 

Figs. 6a, 6b, 6c, and 6d, show the various models employed, and Figs. 

show representative fringe patterns. 


Sources Error the Photoelastic Method 

The results obtained are subject the usual errors the photoelastic 
method, addition the errors peculiar this particular investigation 
which have been mentioned. However, due the special precautions, 
noted below, was felt that the usual errors were kept minimum. 

Time-edge effects the models were practically eliminated keeping 
the models bath Halowax oil. addition, since boundary fringes 
were used basis for calculations, there was cause for such edge 
effects serious. 

Machining stresses were kept minimum making very light cuts 
during the last few passes. The unloaded models showed fringes, and 
uniform intensity light when seen the polariscope. 

The fringe patterns obtained show that creep was probably insignificant, 
since the width the bands was keeping with the stress-gradients. 

The fringe patterns for the models tested the oil bath were somewhat 
fuzzy due certain amount diffusion light the oil. This also 
blurred the edges slightly that was necessary very careful 
measuring the fringe position. 


Test Results from the Constant Thickness Plate 


The value for the difference the principal stresses (p-q) each point 
the plate may computed from Eq. (42). dividing the right side 
this equation obtain the stresses due unit line load, the equa- 
tion becomes, 


(46) 


The position the fringes the z-axis was obtained from the photo- 
graphs and was designated percent the model height. Photographs were 
taken with both dark and light backgrounds for each load, that both the 
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integer Jand half order fringes were obtained. For all tests, the fringes 
were discarded which corresponded maximum shear stress greater 
than the fringe elastic limit, assumed 3500 psi. shear. The values 
(p-q) from tests and theory were then plotted for corresponding values 
percent height, shown Fig. 13. 

From Fig. 13, seen that the test values (p-q) are within about 
the theoretical values the region between and percent model 
height. The values are low percent, high percent, and good 
agreement about percent model height. The values are low near 
the top region the plate due the plastic flow under the pins which 
caused redistribution the stresses. For the lower region the plate, 
the stresses are high due the shape the model, and method support. 


Test Results from the Wedge 


use the analogy, the values h(p-q) the wedge were com- 
pared with the value (p-q) obtained from Borowicka’s solution using the 
value 0.268, (obtained from Eq. (37), using 0.366 for the Bakelite). 
integration Borowicka’s solution give the solution due 
infinite line load, the value for the difference the principal stresses was 
found be, 


This corresponds increase 22.5% over the value obtained for the 
constant modulus elasticity. 

the wedge loaded total load pounds, then the values 
h(p-q) obtained correspond the values (p-q) due line load 
pounds/unit length the point-isotropic material. The values for unit 
line load are again given Eq. (46). 

seen Figs. and 10, that the loading pin has dug into the top 
the wedge considerably, that the load longer concentrated 
point. order correct for this, was assumed that the load the top 
surface was distributed parabolically over the bearing width. was 
desired compare the results for several wedges and loads, and thus dif- 
ferent bearing widths, the correction was applied the experimental 
values rather than the theoretical values. obtain correction factor, 
the ratio (p-q) for parabolic distribution load (p-q) for concen- 
trated load the z-axis was calculated, assuming the radial stress distri- 
bution for and (Table 3). Fig. shows the plot this 
ratio function the angle where the angle subtended the 
load the point under consideration. 

The values (p-q) given Eq. (46) were then divided 
the factor read from Fig. 14, and plotted against percent model height 


shown Fig. 15. The theoretical values (p-q) for the constant modulus 

and for are also plotted this figure. 
The test results for the wedge, corrected for spreading the load, are 

lower than the theoretical solution the upper percent the model 


height. Fig. curve was drawn through the experimental values. 
However, curve through the average the plotted points indicates varia- 
tions from theory (curve and percent 10, 15, 20, and 
percent model height, respectively. The scatter results considerably 
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larger than was obtained using the constant thickness models. This due 
the following three factors: (a) the small loads applied the wedge were 
not measured accurately the loading frame, (b) the fringe count 
much lower than the plate, therefore any residual stresses the model 
become more significant, and (c) the lower stress gradient produces 
broader fringes that location their mid-position more subject 
personal errors. 

These results for (p-q) the axis the wedge are not too satisfactory, 
due the rather large correction factor required account for the distri- 
bution the load. However, additional corroboration the analogy 
was obtained comparing the shape the fringes the shape the 
theoretical fringe. This was done for fringes the region the wedge 
which corresponds the region the constant thickness plate for which 
the fringes are circular. Since the stress distribution for and 
=1/3 not much different than that for 0.268, the fringes the 
wedge were compared with the theoretical shape for For this 
solution, shown Table that principal stress not equal 
zero and result, the expression for (p-q) given as, 


Therefore the theoretical fringes are lines denoting, 
(49) 


Fig. shows theoretical fringe determined Eq. (49) and circle 
which corresponds the theoretical fringe for constant thickness plate. 
The plotted indicate the shape fringes wedge 1.50 0.90 
inches, 0.014 inches thick the top, obtained for loads 20.6 and 41.2 
The experimental and theoretical shapes are seen close 
agreement except very near the top surface, where the spreading the 
load again caused fringe distortion. This correspondence between the 
theoretical and test fringe shapes thus considered acceptable 
confirmation the mathematical anaiogy. 


Conclusions 


The problem the stress distribution semi-infinite body due sur- 
face loading has been analyzed for certain non-homogeneous materials. 
was assumed that the material isotropic any particular point and that 
Poisson’s ratio constant throughout the region, but that the modulus 
elasticity function depth beneath the top surface, although constant 
any given depth. This called point-isotropic material. 

The general equations elasticity for such material have been derived 
for general variation the modulus, and the known solutions due 
Frohlich, Ohde, and Borowicka have been examined with regard their 
validity and limitations. was shown that superposition the correct 
solutions permissible. 

There shown mathematical analogy between the stress distri- 
bution semi-infinite body point-isotropic material due line 
strip loads, and the stress distribution due similar loads semi- 
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infinite plate homogeneous isotropic material. The thickness the 
plate must vary with depth manner corresponding the modulus vari- 
ation the point-isotropic material. Thus becomes possible use the 
photoelastic method determine stresses point-isotropic materials. 

The mathematical analogy confirmed experimentally use 
wedge-shaped photoelastic model. This was analogous point-isotropic 
material for which the modulus elasticity increased linearly with depth 
from zero value the surface. 

Since the variation modulus may cause stress distributions appreci- 
ably different from those obtained for the homogeneous isotropic materials, 
its effect may interpreted concentration.’’ Boundary con- 
ditions, such rigid layers, hydrostatic pressures, also affect the 
stress distributions the overlying material, that necessary 
identify the contribution each the experimental work. was also 
shown that the photoelastic method may applied certain problems 
layered systems provided the results are recognized inaccurate the 
interfaces the layers. 

suggested that the photoelastic analogy may used two ways. 
the first case, complete stress distribution the model may ob- 
tained and compared with the complete stress distribution plate 
similar size. using such comparison, considerations must made 
with regard the method support and the stress concentrating effects 
induced thereby. For the second case, the ratio (p-q) along the z-axis 
the model the value (p-q) along the z-axis similarly supported 
plate determined the region between and percent the model 
height. This ratio may then used obtain appropriate values the 
stress concentration factor, This value may then used 
Frohlich’s equations give approximation for the stress distribution 
the point-isotropic material. The validity this method depends upon 
the relation between and the elastic constants the material. However, 
the approximation reasonable, this method convenient, since 
equations may integrated easily. 

limitation the method that solutions can obtained only for 
values Poisson’s ratio between about 0.23 and 0.33, since presently 
photoelastic materials have Poisson’s ratios varying from about 
0.30 0.50. 
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TABLE 


Vertical Stress Along Loading Axis ina 
Material for which 


Poisson's Ratio 


TABLE 


Vertical Stress, due Point Load Acting Vertically 
the Surface the Semi-Infinite Body 


R* “TT 


Infinite Series Infinite Series 


Infinite Series Infinite Series 
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1/2 1.00 
1/3 1.33 
1.76 
4 
2-2 


TABLE 


Vertical Stress, due Line Load q/unit Length 
Acting Vertically the Surface the Semi-Infinite Body 


*whe — and f 0.7218 
(1/2) 


cos* 
211-25 


Fig. Coordinate Systems 
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(a) Cartesian Coordinates (b) Cylindrical 
three dimensions) (axially symmetric) 


Spherical Coordinates (d) Cartesian Coordinates 
(radial stress distribution) (two dimensions) 


(e) Plane Polar Coordinates 


Fig. Conventions for Positive Stresses 
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Cardboard 


Plexiglas 


Steel Baseplate 


Fig. Method Loading Models 
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Fig. 5. Loading Frame with Model and Tank in Place 
Pp 


Plate, Analogous Semi- (b) Wedge, Analogous Semi- 


Analog Layer Base (d) Analog for Two Elastic Layerson 
Slippery Rigid Base 


Fig. 6. Shapes of Models Tested 
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Fig. 10. Bakelite Wedge (Fig. 6b) Under Load of 72.1 Ibs. 


: 
Fig. Bakelite Wedge (Fig. 6b) Under Load 41.2 
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Fig. 13. Experimental and Theoretical Values of (p - q) along the 
Axis of Loading for the Plate Model 


Fig. 14. Correction Factor for Parabolic Distributions of Load 
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Theoretical Value, 
Top 


Percent Model Height 


Fig. 15. Experimental and Theoretical Values of (p - q) along the 
Axis of Loading for the Wedge Models. (Analogous to 
E = kz for uw = 0.268) 


Load Ibs 
Theoretical Fringe 


Theoretical Fringe 
constant 


Fig. 16. Theoretical and Experimental Fringe Shapes for the Wedge Model 
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